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Abstract 

A method for constructing interferograms, schlieren, and shadowgraphs from ideal- and 
real-gas, two- and three-dimensional computed flow fields is described. The computational grids 
can be structured or unstructured, and multiple grids are an option. The constructed images are 
compared to experimental images for several types of flow, including a ramp, a blunt-body, a 
nozzle, and a reacting flow. The constructed images simulate the features observed in the exper- 
imental images. They are sensitive to errors in the flow-field solutions and can be used to identify 
solution errors. In addition, techniques for obtaining phase shifts from experimental finite-fringe 
interferograms and for removing experimentally induced phase-shift errors are discussed. Both 
the constructed images and calculated phase shifts can be used for validation of computational 
fluid dynamics (CFD) codes. 


Introduction 

For decades experimental interferograms, schlieren, and shadowgraphs have been used 
for quantitative and qualitative flow-field studies. These three images are created by passing 
light through the flow field, and the recorded intensity patterns are functions of the phase shift and 
angular deflection of the light. 

In infinite- and finite-fringe interferograms, the recorded intensity patterns (fringes) are 
caused by phase shifts (optical pathlength differences). These phase shifts result from variations 
in the flow-field density and are proportional to path integrals of the refractive index. The path of 
integration is the path that the light follows through the flow field. 

For two-dimensional (2-D) and axisymmetric flow fields, point information can be 
extracted from interferograms and compared to computational results. The first step in extracting 
this information is to calculate the phase shifts from the interferogram’s fringe patterns. These 
phase shifts can be obtained from either infinite- or finite-fringe interferograms. However, in 
flow-field regions where there are only small changes in the density and, hence, fractions of 
fringe shifts, calculating the phase shifts from finite-fringe interferograms will give more 
accurate results. 

Several methods for calculating the phase shifts from finite-fringe interferograms exist. 
These methods involve either tracing individual fringes, 1 - 2 fitting sinusoidal functions to 
intensity variations in a single interferogram, 2 or using multiple interferograms and phase step- 
ping techniques. 4 After the phase shifts are found, a transformation is applied, and the refractive 
index (and density for non-reacting flows) at every point in the flow field is obtained. This point 
information can then be compared to flow-field solutions. 

For three-dimensional (3-D) flows, the transformation from the experimental phase shifts 
to point information cannot be made, and other methods of comparison must be used. For example, 
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Strike 5 presented theoretical interferograms over 15 years ago, Watkins 6 used path-averaged den- 
sity contours to compare flow-field solutions with experimental infinite-fringe interferograms, 
and Havener 7 tracked individual fringes through flow-field solutions and compared these results 
to experimental finite-fringe interferograms. These methods do allow for comparison of computed 
and experimental flow fields but not necessarily on a one-to-one basis. A graphics package that 
was developed by Tamura 8 and is available in Japan does have the capability of constructing 
interferograms, as well as schlieren and shadowgraphs, from computed flow fields. With this 
package, one-to-one comparisons of computation and experiment are possible. 

In schlieren and shadowgraphs, the intensity patterns (dark and light regions) are 
governed by the angular deflection of the light as it passes through the flow field. The angular 
deflection is proportional to the path integral of the refractive-index gradient. Flow-field features 
which have large refractive-index (density) gradients, such as shocks, shear layers, and expan- 
sion fans, are recorded in these images. 

Both schlieren and shadowgraphs are used for flow-visualization studies and for locating 
flow-field structures. Schlieren contain sharp details of the flow-field structures as well as 
information concerning the direction of the light when it exits the flow. Shadowgraphs also con- 
tain details of the flow field. However, they contain no information as to the direction of the light 
and are generally less sharp than schlieren. When these images are compared to computed flow 
fields, contour plots for only one plane of data (usually the symmetry plane) are typically used. 
The choice of contour levels can mask some of the features observed in the experimental images. 
In addition, if the flow is 3D or the model is free to roll, no single computational plane provides all 
the information necessary for realistic comparisons to experimental images. 

The information required for constructing interferograms, schlieren, and shadowgraphs 
is contained in the flow-field solutions. The refractive index is proportional to the gas-species 
densities, and these densities are always calculated in the solutions. By developing routines to 
evaluate and integrate appropriate functions of the refractive index, direct comparisons of 
computed and experimental results can be made. 

In this paper, the CISS (Constructed Interferograms, Schlieren, and Shadowgraphs) code is 
described. CISS is software that constructs images from ideal- and real-gas, 2-D, axisymmetric, 
and 3-D computed flow-field solutions. The computational grids can be structured or unstructured, 
and multiple grids are an option. This software differs from Tamura s 8 in that with CISS the 
computational cells are not divided into tetrahedra; CPU time requirements are reduced by 
eliminating the need for ray tracing; all image points are calculated and no interpolation between 
image points is required; the integration arguments are functions of the refractive index, not the 
density, and images for real-gas flows can easily be constructed; and the method for constructing 
shadowgraphs is different. Constructed and experimental images are shown for several types of 
flow fields, and the effects of grid resolution and solution convergence are discussed. 

In any comparison of experimental and constructed images, the issue of optical 
aberrations should be addressed. These aberrations must be determined individually for each 
experiment, and they should be removed from the experimental image before comparisons between 
experiment and computation are made. Although a general treatment of this subject is beyond the 
scope of this paper, an experimental finite-fringe interferogram is used to demonstrate the 
importance of identifying and removing aberrations. In this example, a method for extracting 
phase shift information from the experimental interferogram is described. The phase-shift errors 
that are caused by optical aberrations are then identified and removed from the calculated phase 
shifts. 
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Integrating these functions along a line of sight gives the total angular deflections, £x and £y. 

In shadowgraphs, no knife edge is used and the dark and light regions are caused by the 
concentration and divergence of light (Fig. 2). In many texts (e.g. Refs. 9 and 10), the function 
used to describe shadowgraphs is the line integral of the second derivative of the refractive index. 
In CISS, the function of the refractive index used for constructing shadowgraphs has been chosen to 
be equal to that used for schlieren. At the exit of the computed flow field, the angular deflection is 
known. If the distance to the image plane is also known, then the area on the image plane where 
the deflected light falls is known, and the shadowgraph can be constructed by adding the contribu- 
tion of the deflected light beams at each point on the image plane. It should be noted that if 
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there will be no convergence or divergence of the light and, hence, no variation in shade. The form 
of this equation provides the motivation for describing shadowgraphs as line integrals of the sec- 
ond derivative, although this description is not precisely correct. Tamura 8 constructed shadow- 
graphs by calculating the divergence (V ) of the integrated density gradient. In CISS the light is 
actually projected onto the image plane. The advantage of constructing shadowgraphs in this 
manner is that the thickness of the dark and bright regions can be controlled by the placement of 
the image plane. This dependency of the dark and bright regions on the image plane's position is 
observed in experimental shadowgraphs. 

Integration Scheme 

The intensity variations observed in the images are governed by integrals of the 
refractive-index functions described in the preceding section. If an image is divided into an MxN 
array of pixels, these integrals must be evaluated for every pixel in the image. The integral for the 
(m t n )th element of the array is 


F m n = J* / [ft (Xm(Q’yn(0> o] 

0 

where ds is the element of arc length. The light's path, defined by (x m (Q,y n ( Q>Q for 0 < £<z, 
intersects the image plane at the (m,n)th element. This integral can be written as 


N 2l+i 
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where the i summation is over all computational cells which contain segments of the line of 
integration. The interval z^ < f <Zj + j defines the intersection of the line of integration with the i th 
computational cell (Fig. 3). 

Tracing the actual light path as it bends through the flow-field solution and integrating the 
appropriate function of the refractive index along this path is computationally expensive. The 
computer resources required for producing constructed images can be appreciably reduced by 
approximating the light’s path by a straight line perpendicular to the image plane. With a straight 
line approximation, x m and y n become independent of £, and F m n can be rewritten as 
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F m,n = X J f ^ X m,y n ,0\ d£ 

i=l 2 . 

The order of evaluation and summation of the integrals has no effect on F m n and, therefore, no 
ray tracing is required. It is necessary only to find the intersection of the line of integration with 
each computational cell, calculate the integral for this segment, and add it to the appropriate sum. 

To determine whether the line of integration intersects the computational cell, it is 
necessary to determine if the point ix m ,y n , £) falls within the computational cell for any value of C- 
In CISS, the test procedure is as follows. Each computational cell (for example, the cell shaded in 
gray in Fig. 3) is defined by eight points,^* = tei,y„ Zj), and six surfaces. Each surface is then 
described by two triangles, and for each triangle, three cross products are formed: 

(£-£»> x ] 

(£-3m > X > 

In these equations, = (x m j nr z) gives the position of the (m,n)th pixel of the image array, and 
the^ = (x;,y;, z) for i = a, b and c are the projections onto the image plane of the vectors that define 
the vertices of the triangle. To minimize CPU time, testing is performed only for x m and y n that lie 
within the minimum and maximum x and y values for the computational cell. The sign for all 
three cross products is the same if and only if the point lies within the triangle, and the triangular 
surfaces through which the line of integration enters and exits the computational cell can be deter- 
mined. Once these surfaces are determined, interpolated values for Zj, z, + j, f(x m , y n , zj, and 
f(x m , y n > Zi+i) are obtained, and the line integral of the function / is evaluated using the 
trapezoidal rule. In CISS, three functions of the refractive index are integrated simultaneously: n- 
n 0 , (l/n)dnl dx, and (l/n)dn/dy. 

Combining this searching routine with the straight line approximation minimizes the 
CPU time required for creating the images. For regions of the flow without shocks, the straight line 
approximation should have a minimal effect on the constructed images. In shock regions, the re- 
fractive index changes rapidly, and the approximation may introduce errors; however, these 
errors are in many cases on the order of, or smaller than, the solution and experimental errors. 

Estimation of Errors 

The constructed images have several sources of errors. These include errors in the flow- 
field solutions, the approximation of the light path by a straight line, the interpolation schemes, 
and the accuracy of the integration method. Although solution errors are of great importance, they 
will be mentioned in this paper only when they directly affect the process of constructing images 
from flow-field solutions. 

The error, E, introduced by the straight line approximation can be written as 

E =j f(x(Qy(Q,Qds - J7(O,O,0dC = \f (x(CW0»0M» * <*C> + JL f (*(OM.O ■ fW,0,Q]dC 

For simplicity, the x,y-coordinates of the image pixel have been set equal to zero. The first integral 
on the right hand side of this equation represents the error introduced by the difference in path 
length; the second term represents the error introduced by approximating / by its value on the 
straight line. By writing the element ds as 
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and by expanding the integrands in terms of /, x, y, dx/dz, and dy/dz, the expression for the error 
can be written as 


E = = 
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If the light and integration paths initially coincide, x ~ jdx/dz d£ and y - \dy/dz d£. The 

derivatives dx/dz and dy/dz can be approximated by + JV dp! dx dC, and £y 0 + jKdp/ dyd( where 
e Xo and £y 0 are the initial angles between the light and integration path. Substituting these 
expressions into the above equation gives 
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where C is the length through the disturbed region of the flow field. 

For interferograms, the function f controls the phase shift and is given by / = 2n(n - n 0 )/A = 
2m r(p - p 0 )/ A; A is the wavelength of the light used to create the interferogram. The estimated error 
in the calculated phase shift reduces to 
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Under most conditions, k p is much less than one, and the above estimate can be manipulated to 
give 0(A<p), the order of the error in the calculated phase shift: 


0(A<p) ~ O 




In schlieren and shadowgraphs, the light and dark regions are governed by the deflection 
of the light through the flow field. The function that controls this deflection is proportional to the 
gradient normal to the light path of the refractive index (and hence the density). The order of the 
error in the calculated deflection angle £ x is 
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A similar equation can be obtained for £y. 

For all three images, the errors introduced by the straight line approximation are functions 
of the initial angle between the integration path and the actual light path. This initial angle, for 
instance, might be caused by the deflection of the light as it passes through an outer shock. To 
estimate this angle, consider light passing through a 2-D shock. The deflection angle e Q is given 
by 


e n = sin * 


cos0 


\l 1 - (— 'f cos 2 0 - — sin0 
y l ^2 J n 2 


P 


The angles 9 and £ 0 are defined in Fig. 4; and rt 2 are the refractive indices on the incident and 
transmitted sides of the shock. The maximum deflection angle through the shock is 


£o = sin 


n\ - 

nl 


+ nj J 


k(A p ) 


where A p is the change in density through the shock, and the initial angles in the error estimates 
can be replaced by k Ap . 

The estimates on the order of the errors introduced by the straight line approximation are 
functions of the flow-field density, its derivative, and the Gladstone-Dale constant for the gas. For 
air, Kp s tp is approximately 3xl0‘ 4 where p stp is the density at standard room temperature and 
pressure. Normalizing the density p by p s t p allows the error estimates to be rewritten as 


and 
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In these expressions and A have been nondimensionalized by t 

The wavelength of light used to create most interferograms is on the order of 4-7xl0‘ 7 m, 
and t! A can be a large number. If f is on the order of meters, the order of the error in the calculated 
phase shift will be 


0(A<p) - 0 
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For the straight line approximation to be valid for flows of a few meters, the density jumps must be 
less than 10 times the density at standard room temperature and pressure and the normalized 
density gradients in regions of the flow away from shocks must be less than one. It should be noted 
that if these conditions are violated, areas of the interferogram can still be accurate. These areas 
include regions where the light is nearly normal to the shock and where the light does not pass 
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through regions of rapidly varying density. The schlieren and shadowgraphs will be accurate for 
a much wider range of conditions. 

Two other sources of error are caused by inaccuracies in the interpolation and integration 
methods. As described earlier, the surfaces of the computational cell are triangulated, and the 
triangle through which the light beam enters or exits the computational cell is determined. A 
linear interpolation then uses values at the three vertices of the triangle to calculate the z- 
coordinate, the refractive index, and the refractive-index gradient at the entrance or exit point. 
Once these values are calculated, the integral is found using the trapezoidal rule, a second-order 
method. Both the linear interpolation and the integration method can introduce errors. However, 
these errors are minimized when the grid has been adapted to the flow-field solution and large 
gradient regions are well resolved by the computational grid. 

Construction of the Images 

For infinite-fringe interferograms, the intensity pattern is simply proportional to the sine 
of the phase shift <p: 



j (n - n 0 ) dC 


For finite-fringe interferograms, a linearly-varying phase is added: 


5<p = k x x + k y y 


where k x and k y are functions of the angle between the reference and object beams. 

For schlieren, the intensity is a function of the angular deflection. As an example, 
consider a horizontal knife edge (Fig. 1) and a rectangular section of the initial light source with 
height h. As this rectangular section passes through the flow field, it is deflected and exits the flow 
field at an angle. When this light is focused at the plane of the knife edge, the rectangular element 
has a height of h' and has been deflected vertically by a distance yey] y is a function of the distance 
between the flow field and the knife edge. There is also a horizontal deflection; however, it is 
unimportant in this case. If the horizontal knife edge is placed at the center of the plane, the ratio of 
light not blocked by the knife edge to the amount of the original light is 
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and the intensity patterns are functions of this ratio. 

As mentioned previously, shadowgraphs are not constructed in CISS from line integrals of 
second derivatives of the refractive index. Instead they are given by variations in intensity 
caused by the deflection of the light. In CISS it is assumed that the distance through the region of 
disturbed flow is much less than the distance to the image plane, and that whereas the deviation of 
the light path from a straight line path is not significant in this region, it is large enough to cause 
variations in intensity at the image plane. The shadowgraphs are constructed by tracking 
elements of light after they exit the disturbed flow region. The location and size of each light 
element is governed by four light rays which define the corners of that element, and the direction 
of these four light rays is given by the angle at which they exit the disturbed flow region, that is, 
shadowgraphs are also governed by £ x and £y . Using these angles, the size and position of each 
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light element at the image plane is calculated, and the amplitude of the electric field, E mny at the 
(m,n)th pixel in the image plane is increased by 


AEmn = Eo 


Cf r)A r 
Ct 


Here, E 0 is the amplitude of the initial field, Cf r\A mn is the area of the intersection of the deflected 
element with the (m,n )th pixel of the image plane, and C; is the initial area of the deflected 
element. The intensity pattern observed in the shadowgraph is proportional to the square of E mn . 

The time requirements for constructing these images from flow-field solutions are 
linearly proportional to the grid size and the number of pixels defining the image. For 2-D flow 
fields, the image construction takes only a few seconds on a CRAY-YMP. For 3-D and 
axisymmetric flow-field solutions, more CPU time is required. For example, the construction of 
an image with 400,000 pixels from a flow-field solution with 250,000 grid points will take a few CPU 
minutes. 


Examples of Constructed Images 

Representative images are shown in this section for several types of hypersonic flows 
including a nozzle, a blunt-body, a ramp, and a reacting flow. The computed flow fields include 2- 
D, axisymmetric, and 3-D solutions. The purpose here is not to provide a detailed analysis of the 
computations but to simply demonstrate the ease of using constructed images for comparing 
experimental and computed results. 

In Fig. 5, experimental and constructed infinite-fringe interferograms are shown for a 2- 
D nozzle flow. 11 The experiment was performed in the Electric Arc Shock Tunnel at NASA Ames 
Research Center, and the test gas was nitrogen. The test conditions were such that nitrogen disso- 
ciation and recombination were important, and a real-gas CFD code was used for the flow-field 
solutions. Note that by constructing an infinite-fringe interferogram from the flow-field solution, 
one-to-one comparisons of computation and experiment are possible. The areas along the center 
line where the experimental and computed densities agree or disagree easily can be seen by the 
matching (on the left hand side of the interferogram) or mismatching (on the right) of the fringes. 
Differences between the computed and experimental boundary layer thickness and the density 
near the wall are readily seen from the differences in the shapes and orientations of the in- 
dividual fringes. 

In Fig. 6, constructed and experimental infinite-fringe interferograms are shown for a 
ballistic range test of the Aero-assist Flight Experiment (AFE). The flow-field solution was 
calculated using an ideal-gas CFD solver; 12 the grid was adapted to the flow-field solution using 
SAGE. 13 A one-to-one comparison of the two interferograms is not possible; the flow-field solution 
is for a model with no yaw, the experimental model did yaw. However, there are similarities. The 
constructed and experimental interferograms both show similar light and dark patterns, and the 
number of fringe shifts is comparable. The location of the fringes provides information as to the 
values of the computed and experimental integrated refractive indices. 

In Fig. 7, constructed and experimental schlieren are shown for a 2-D, 15° compression 
corner, the flow is from left to right. The flow-field solution was computed using an ideal-gas 
solver. 14 For both schlieren, a horizontal knife edge was used. The shocks, shear layers, and sepa- 
ration region observed in the experimental schlieren are simulated by the constructed schlieren. 
The dark regions in the constructed schlieren are caused by rapid increases in density across the 
computed shocks, and they correspond to dark regions in the experimental schlieren. The light 
regions in both the constructed and experimental schlieren are identified with boundary layers, 
and they are caused by density gradients near the surface. A shear layer (identified with a light 
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region) lifts off the surface when separation occurs, and the extent of the computed and experi- 
mental separation regions can be compared. The grid for the computational study was not adapted 
to the flow-field solution. The sharpness of the features in the constructed schlieren could be fur- 
ther improved by using grid adaptation to reduce solution errors. 

In Fig. 8, constructed and experimental shadowgraphs for a generic, National Aero-Space 
Plane (NASP) nozzle configuration, the Single Expansion Ramp Nozzle (SERN), are shown. Both 
the freestream and nozzle flows are supersonic. The flow-field solution in Fig. 8a is a 2-D, ideal- 
gas solution at the symmetry plane; the solution in Fig. 8b is 3D. 16 The interaction of the external 
flow with the nozzle plume is complicated, and all of the features observed in the experimental 
shadowgraph are captured in the constructed shadowgraphs. These features include the leading 
edge shock, separation from the cowl, and the shocks, shear layers, and expansion fans caused by 
the interaction of the nozzle plume with the freestream. As with the schlieren, dark and light 
regions in the constructed shadowgraph correspond to dark and light regions in the experimental 
shadowgraph. However, whereas the shocks in the schlieren are described by either a single dark 
(or light) line, shocks in the shadowgraphs are described by a combination of dark and light lines. 
It should be noted that there are two parameters in CISS that allow the effective sensitivity and 
exposure of the film to be modeled. In this case, these parameters were chosen to emphasize the 
location of the flow-field features; they do not simulate the experimental film characteristics. 

In the two constructed shadowgraphs in Fig. 8, the effects of grid resolution and solution 
errors can be seen. The grid in the 2-D solution is very fine, and the constructed shadowgraph 
resembles the experimental shadowgraph. For the 3-D solution, although multiple grids were used, 
the gridding was much coarser than that used for the 2-D solutions. The shocks, shear layers, and 
expansion fans were not as well resolved in the 3-D flow-field solution and, hence, are more 
diffuse in the constructed shadowgraph. The leading edge shock in the 3-D computation has 
become so diffuse that it is barely detectable in the constructed shadowgraph. 

In Fig. 9, a constructed shadowgraph is shown for an axisymmetric, real-gas, flow-field 
solution of a ballistic range shot of a hemisphere cylinder into a combustible mixture . 11 • 18 Also 
shown is an experimental shadowgraph at slightly different test conditions. 19 This flow is very 
complex, and many of the features observed in the shadowgraph are 3D and difficult to identify in 
ordinary contour plots. In the constructed shadowgraph most, if not all, of the features in the exper- 
imental shadowgraph are captured. It should be noted that in most shadowgraphs, the leading edge 
of a shock is usually seen as a dark line immediately followed by a bright line. In this case, the 
optical set up of the experiment reversed the order of the dark and bright lines. CISS was able to 
duplicate this effect by using a negative number to represent the distance to the image plane, and a 
one-to-one comparison of experiment and computation is possible. 


Interpreting Finite-Fringe Interferograms 

In the previous section, it has been demonstrated that the ability to construct images from 
computed flow fields facilitates the comparison of experiment and computation; hence, the CISS 
code provides a useful tool for CFD code validation. The capability for removing optical 
aberrations and for extracting fringe-shift information from finite-fringe interferograms is also 
invaluable for code validation. The interference patterns recorded in interferograms are 
functions of the phase shifts caused by variations in the refractive index. For low density flows, 
there may be only fractions of fringe shifts in the regions of interest (for instance, the wake 
region), and the phase shifts cannot be determined from the fringe patterns in infinite-fringe 
interferograms. However, they can be determined from the patterns in finite-fringe interfero- 

In finite-fringe interferograms, a linearly-varying phase is added to the phase shifts 
caused by the refractive-index variations. This results in straight, parallel, and equally-spaced 
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fringes in the freestream regions of the interferogram. In regions where the flow is disturbed, the 
fringes bend, converge, and/or diverge. 

In the past, a method for determining the phase shifts involved tracing individual fringes. 
By measuring the displacement of individual fringes from the freestream pattern, the phase shift 
could be determined. The resolution of this method is limited by the fringe spacing, and it is time 
consuming. 

By digitizing the finite-fringe interferogram, gray-scale information (the intensity 
pattern) can be used to calculate the phase shift at every point in the interferogram. Different meth- 
ods for finding the phase shift from digitized interferograms have been developed, some of them 
requiring multiple finite-fringe interferograms 4 and others requiring a single interferogram. 3 
Included here is a description of a technique that uses a single digitized interferogram and that 
fits sine waves to gray-scale values. 

When calculating the phase shift from finite-fringe interferograms, the frequency and 
phase of the freestream fringes first must be determined. A region of the digitized interferogram 
is identified with the freestream, and the freestream intensity pattern is approximated by 

I 0 = A 0 + B 0 sin( iqjc + <t>o) 

The freestream spacing, which is determined by k xo , and the phase, <po, are found by an iterative, 
least-squares procedure. 

Once k xq and <p 0 are known, the phase shifts are determined pixel by pixel. First, a single 
pixel is identified in the digitized interferogram and the gray-scale values are obtained along a 
line that is centered about this pixel and that is perpendicular to the freestream fringes. The length 
of this line is equal to one freestream wavelength. The intensity pattern is then approximated by a 
sinusoidal function 


I = A + B sinfxy: + <p) 

and one of three methods is used to calculate the four parameters in this equation. In the first 
method, is assumed to be approximately equal to x^ 0 , and three points that are distributed over the 

length of the line are used to calculate A, B, and <p. In the second method, the frequency is not 
assumed to be equal to the freestream frequency, and four points that are distributed over the length 
of the line are used to calculate A, B, x^. , and <j>. In the third method, all of the points along the line 
are used, and A, B, k x , and <p are found using an iterative, least-squares fitting procedure. Once 
k xq and <p are known, the phase shift, 


5<j> = (K x - x^ 0 )x + (<t>- 0 O ) 

is calculated. A new pixel in the interferogram is then chosen, and the process is repeated until the 
phase shift is determined for every pixel in the digitized finite-fringe interferogram. The differ- 
ences in the calculated phase shifts using these three methods are minimal, and the first two 
methods are computationally much faster. 

In Fig. 10a an experimental finite-fringe interferogram is shown for a blunt cone. This 
interferogram was taken in the ballistic range at NASA Ames Research Center. The test section 
pressure was 0.08 atm. 20 In Fig. 10b, the calculated phase shift modulo 2n is plotted. It should be 
noted that jumps from white to black represent only small changes in phase. The freestream phase 
shift in the majority of the lower half of Fig. 10b is white or black; hence, the calculated phase is 
nearly constant in this region. However, in the lower right hand corner and in the upper half of the 
plot, there are variations in the calculated freestream phase shift. Investigation of the digitized 
image did indicate that the fringe spacing was actually varying in the freestream. Since the 
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freestream flow in the ballistic range was uniform and quiescent, these variations are errors that 
have been introduced by the experimental apparatus. 

The finite-fringe interferogram in Fig. 10a was created using a dual plate, holographic 
interferometer. 20 Variations in the freestream fringe spacing could have several sources 
including the spatial quality of the laser, the optics, and the reconstruction process. In order to 
remove these variations, a least-squares procedure was used to fit polynomial functions to the 
freestream phase. These functions were then used to correct the phase shift for the entire image, 
and the results are shown in Fig. 10c. In this plot, the freestream region is nearly all white and 
black, and the variations in the calculated freestream phase shifts are now very small. The esti- 
mated errors in the phase shift calculations are less than 10%. 

Differences in the wake regions for the uncorrected and corrected calculated phase shifts 
should be noted. The uncorrected phase-shift calculations indicate that there is a neck region in 
the wake flow. Once the freestream phase variations have been removed, the neck region disap- 
pears. Qualitative as well as quantitative analyses of finite-fringe interferograms require re- 
moval of experimental errors such as those caused by optical aberrations. 


Concluding Remarks 

The CISS code has proven its capability in constructing interferograms, schlieren, and 
shadowgraphs from a variety of ideal- and real-gas, 2-D, axisymmetric, and 3-D flow-field 
solutions. For the examples shown here, the constructed images are very similar to the 
experimental images - the flow features in the constructed images are represented by the same 
variations in light and dark as in the experimental images. The fringe patterns in the constructed 
finite- and infinite-fringe interferograms simulate the patterns in the experimental interfer- 
ograms, and one-to-one comparisons of constructed and experimental interferograms can be 
made. The shocks and shear layers in the constructed schlieren and shadowgraphs are defined by 
the same intensity variations as in their experimental counterparts. In addition, 3-D effects are 
included in the construction. Hence, one-to-one comparisons of constructed and experimental 
schlieren and shadowgraphs also can be made. 

The quality of the flow-field solution does affect the quality of the constructed images. 
Insufficient or improperly placed grid points can cause spreading of shocks and shear layers in 
the constructed images. CISS not only provides an excellent tool for comparing flow-field solu- 
tions with experiment, it can be used to identify problem areas in the flow-field solutions. 

By constructing images from flow-field solutions and by comparing these images to 
experimental images, the ability of the flow solver to simulate experimental flows can be assessed. 
More detailed information for code validation can be obtained by determining phase shifts from 
experimental finite-fringe interferograms and by comparing these phase shifts to computed phase 
shifts. Methods for extracting phase shifts from finite-fringe interferograms do exist. It is 
imperative that the experimental errors caused by optical aberrations be removed from the phase 
shift calculations before either qualitative or quantitative comparisons are made. The ability to 
compare computations and experiments on a one-to-one basis and to extract detailed information 
from finite-fringe interferograms is invaluable for CFD code validation. 
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Collimated Image 



FlowfieW Lens Knife Lens lntenslt >' 
edge pattern 

Fig. 1 Schematic for constructing schlieren. 



Fig. 2 Schematic for constructing shadowgraphs. 




Fig. 4 Refraction of light through a shock. 



Fig. 5 Experimental and constructed infinite- 
fringe interferograms for 2-D nozzle 
flow. 11 






Fig. 6 



Constructed and experimental infinite-fringe interferograms for an AFE model. Three- 
dimensional flow field solution: Venkatapathy; 12 experiment performed at NASA Ames 


Research Center's ballistic range (unpublished data). 



Constructed schlieren Experimental schlieren 

Fig. 7 Constructed and experimental schlieren for a two-dimensional compression ramp. Flow 
solution: Tokarcik et al.; 14 experiment: Delery and Coet. 15 



Fig. 8 


a) 2-D solution b) 3-D solution c) Experiment 

Constructed and experimental shadowgraphs for SERN. Two-dimensional flow solution and 
ixperiment: Ruffin et al. 16 









Constructed shadowgraph Experimental shadowgraph 

Fig. 9 Constructed and experimental shadowgraphs for a hemisphere-cylinder in a combustible 
mixture. Flow field solution: G. Wilson; 1 ^’ 1 ** experiment: Lehr. 19 




Fig. 10 Finite-fringe interferogram for a blunt cone; 

a) experimental; b) phase shift modulo 2ir, c) phase shift 
modulo 2n with freestream variation removed, d) color bar. 


